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Abstract 

Exact Heisenberg operator solutions for independent 'sinusoidal coordinates' as 
many as the degree of freedom are derived for typical exactly solvable multi-particle 
quantum mechanical systems, the Calogero systems based on any root system. These 
Heisenberg operator solutions also present the explicit forms of the annihilation- creation 
operators for various quanta in the interacting multi-particle systems. At the same time 
they can be interpreted as multi-variable generalisation of the three term recursion re- 
lations for multi-variable orthogonal polynomials constituting the eigenfunctions. 

1 Introduction 

Modern quantum physics is virtually unthinkable without annihilation-creation operators, 
which are defined as the positive/negative energy parts of the Heisenberg field operator 
solutions of a free field theory, an infinite collection of independent harmonic oscillators. 
These annihilation-creation operators map an eigenstate of a free Hamiltonian into another, 
but not connecting those of a full theory. Our knowledge of the Heisenberg operator solutions 
of a full interacting theory, on the other hand, is quite limited in spite of the central role 
played by the Heisenberg operator solutions in field theory in general. In the so-called exactly 



solvable quantum field theories, factorised S'-matrices and some of the correlation (Green's) 
functions are the highest achieved points up to now. 

Following our embryonic work pQ on the construction of exact Heisenberg operator solu- 
tions for various degree one quantum mechanics, we present in this paper a modest first step 
in the quest of deriving exact Heisenberg operator solutions for a family of interacting multi- 
particle dynamics. These are the Calogero systems [2], which are integrable multi-particle 
dynamics based on root systems. For the theories based on the classical root systems, the 
A, B, C and D series, the number of particles can be as large as wanted, but not infinite 
as in field theories. A complete set of exact Heisenberg solutions for 'sinusoidal coordinates' 
[U, [3] , as many as the degree of freedom, is derived elementarily in terms of the universal 
Lax pair [U [5j [6] , which is a well-established solution mechanism for classical and quantum 
Calogero-Sutherland-Moser systems [2J, [8] based on any root system [9] . Explicit forms 
of various annihilation-creation operators are obtained as the positive/negative energy parts 
of the Heisenberg operator solutions. They map an eigenvector of the full Hamiltonian into 
another. These sinusoidal coordinates and the corresponding annihilation-creation operators 
provide multi-variable generalisation of the three term recursion relations [I] of orthogonal 
polynomials constituting the eigenfunctions of the Calogero Hamiltonian. 

This paper is organised as follows. In section two, the rudimentary facts and notation 
of Calogero systems based on any root system are recapitulated together with the universal 
Lax pair matrices. In section three the complete set of exact Heisenberg operator solutions 
is derived quite elementarily based on generating functions constructed from the universal 
Lax matrices. Remarks on the special features of the D-type theories are given at the end 
of the section. The final section is for a summary and comments for further research. The 
Appendix gives the list of the preferred sets of weight vectors for explicit representations of 
Lax matrices based on the exceptional and the non-crystallographic root systems. 

2 Calogero Systems 

In this paper we will derive exact Heisenberg operator solutions for the Calogero systems. 
They are one-dimensional multi-particle dynamics with inverse (distance) 2 potential inside a 
harmonic confining potential. They have a remarkable property that they are exactly solvable 
at the classical [10] and quantum [21 [6] levels. The exact quantum solvability has been shown 
in the Schrddinger picture, as the entire energy spectrum is known and the corresponding 
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eigenf unctions can be constructed explicitly by a finite number of algebraic processes (the 
lower triangularity of the Hamiltonian in certain basis) [6]. The exact quantum solvability 
in the Heisenberg picture of the Calogero systems will be demonstrated in the next section 
by constructing the explicit Heisenberg operator solutions for the independent 'sinusoidal 
coordinates' as many as the degree of freedom. Let us briefly recapitulate the essence of the 
quantum Calogero systems [21 [6] together with appropriate notation necessary in this paper. 

2.1 Quantum Hamiltonian 

A Calogero system is a Hamiltonian dynamics associated with a root system [H |9J A of rank 
r, which is a set of vectors in IR' r with its standard inner product. Its dynamical variables are 
the coordinates {qj} and their canonically conjugate momenta {pj}, satisfying the canonical 
commutation relations^ 

fehPfc] = ^jfc, [qj,Qk] = [Pj,Pk] = 0, j,k = l,...,r. (2.1) 
These will be denoted by vectors in W, 

r r r 

def / \ def / -. def \ > o def \ v 2 2 def \ v 2 /n r»\ 

q = (q u ...,q r ), p = (pi, . . . ,p r ), p ■ q = ^PjVji P = Z^Pj> Q = l^ q r ^ 22 ' 

j=i j=i j=i 

The momentum operator pj acts as a differential operator 

• 9 ■ 1 
Pj = , J = 1, . . . , r. 

d qj 

The 'factorised' Hamiltonian is 

«fc.>*&(»-©(«+0 w) 

3=1 J J 

i 2 , u\ 2 | 1 y 9\p\{g\ P \ -i)IpI 2 F (9A , 
~2 P + Y q + 2 \ — G^P °' ( } 

W A (?) = -|g 2 + g\ p] \og\p-q\, g ]p] >0, ^>0. (2.5) 
peA+ 

The summation is over the set of positive roots A+, with A = A + U (— A+). The real 
positive coupling constants g\ p \ are defined on orbits of the corresponding Coxeter group, i.e. 



1 For the A- type theory, it is customary to consider A r _i and to embed all the roots in R' r . This is 
accompanied by the introduction of one more degree of freedom, q r and p r . The genuine A r _\ theory 
corresponds to the relative coordinates and their momenta, and the extra degree of freedom is the center of 
mass coordinate and its momentum. 
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they are identical for roots of the same length. The Hamiltonian is invariant under the finite 
reflection group (Coxeter group, or Weyl group) generated by the set of roots A: 

H(s a (p),s a (q))=H(p,q), Ma G A, (2.6) 

with the reflection s a defined by 

s a [x) =x- (a v -x)a, Vx G W , a y = 2a/\a\ 2 . (2.7) 

Obviously the Heisenberg equations of motion for the coordinates {qj} are trivial 

d 

1>[H, Qj] = j t Qj =Pj, J = 1, • • • , r, (2.8) 
whereas those for the canonical momenta {pj} have the same form as the Newton equations: 

The hard repulsive potential ~ l/(p • q) 2 near the reflection hyperplane H p = f {q G M. r \ p-q = 
0} is insurmountable at the quantum level as well as the classical. Thus the motion is always 
confined within one Weyl chamber. This feature allows us to constrain the configuration 
space to the principal Weyl chamber (II: set of simple roots) 

PW = {q G W\ a-q > 0, a G II}, (2.10) 

without loss of generality. 

The positive semi-definite form of the factorised Hamiltonian (12.31) simply allows the 
determination of the ground state wavefunction: 

HMi) = o, Mi) = z w{q) = II \p ■ i\ 9M ■ e " f ?2 ' ( 2 - n ) 

which is real and obviously square integrable 

(f>l(q)d r q< oo. (2.12) 



'PW 

The constant part E Q of the Hamiltonian (12.31) is usually called the ground state energy 



pGA H 
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The excited energy spectrum is integer spaced and is independent of the coupling constants 
{9\ P \} ®- 



H<p n (q) = £ n 0n(<?), n = f (m, • • • , n r ), Uj eN = Z> , 

r 

£ n = uN n , N n = ^rij fj. 

3=1 

Here {fj} are the integers related to the exponents {ej} of A: 

fj = 1 + ej. 



(2.14) 
(2.15) 



(2.16) 



They indicate the degrees where independent Coxeter invariant polynomials exist. The set 
of integers 

F A = {fuf2,...,fr} (2.17) 
is shown in Table I for each root system A. 



A 


Fa 


A 


Fa 


4-i 


{2,3,...,r,l} 


E 8 


{2,8,12,14,18,20,24,30} 


B r 


{2,4,6,. ..,2r} 


F A 


{2,6,8,12} 


C r 


{2,4,6,. ..,2r} 


G 2 


{2,6} 


D r 


{2,4,...,2r-2,r} 


h(m) 


{2,m} 


E 6 


{2,5,6,8,9,12} 


H 3 


{2,6,10} 


E 7 


{2,6,8,10,12,14,18} 


#4 


{2,12,20,30} 



Table I: The set of integers Fa = {fi, f 2 , ■ ■ ■ , f r } for which independent 
Coxeter invariant polynomials exist E 

Excited states eigenf unctions have the following general structure: 

0n(<?) = Ml) P n(q), 

in which P n {c[) is a Coxeter invariant polynomial in {qj} of degree iV n . 



(2.18) 



2 For A r _i root system, f r — 1 corresponds to the degree of freedom for the center of mass coordinate. 
The B r and C r Calogero systems are equivalent. 
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2.2 Quantum Lax Pair 

The derivation of the Heisenberg operator solutions depends heavily on the universal Lax 
pair which applies to any root system. The universal Lax pair operators [5j [6] are 

L(p, q)=p-H + X(q), X(q) d = f i £ g M ( -^s p , (2.19) 

P eA+ P q 

«M f -^»^«;4EM^x/, (2.20) 

in which / is the identity operator and {s a : a G A} are the reflection operators of the root 
system. They act on a set of W vectors 1Z = f {fi^ G W\k = 1, . . . , d }, permuting them 
under the action of the reflection group. The vectors in 1Z form a basis for the representation 
space V of dimension d. The operator M satisfies the relation [5l [6] 

which is essential for deriving quantum conserved quantities and annihilation-creation op- 
erators. The matrix elements of the operators {s a : a G A} and {Hj : j = 1, . . . ,r} are 
defined as follows: 



def r r / rr \ def 



S^spC) = ^p(m)> (Hj)ija> = Vjbpv, P G A, (i,u eK. (2.22) 



The Lax equation 

i[W,£] = ^L=[L,M] (2.23) 

is equivalent to the Heisenberg equation of motion for {qj} and {pj} for the Hamiltonian 
(12 .4p without the harmonic confining potential, i.e., u — 0. It should be emphasised that the 
l.h.s. of (12.231) is a quantum commutator, whereas the r.h.s. is a matrix commutator as well 
as quantum. The full Heisenberg equations of motion with the harmonic confining potential 
read 

i[H, L*] = j t L± = [L*, M] ± iwL*, (2.24) 
in which M is the same as before (12.201) . and L ± and Q are defined by 

L ±d =L±iujQ, Q = q-H, (L+) f = L~, (2.25) 
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with L, H as earlier (12.191) . (I2.22p . One direct consequence of the Lax pair equation (12.241) 
is the existence of a wide variety of quantum conserved quantities: 

d k 

— Ts(L £l L €2 ■■■L £k ) = 0, Ej G {+,-}, V/eG2N, ^e i = 0, (2.26) 

Lib 

in which Ts(A) denotes the total sum of a matrix A [5, 6J with suffices /j,,v ETZ: 

Ts(A) ^ ]T A, u . (2.27) 

This is a simple outcome of the property (12.211) of the matrix M [H EJ [HI [12]. The quantum 
conserved quantities are the simplest example of more general results: 

i[H, Ts(L ei L e ' 2 ■ ■ ■ L ek )} = jTs(L ei U 2 ■■■L ek ), Vk G N, 

k 

= iumTs(L ei L e2 ---L^), ^e i = m. (2.28) 

i=i 

In other words, Ts(L ei L e2 ■ ■ ■ L €k ) shifts the eigenvalue of 7i by m unit (or, um) when it acts 
on any eigenstate of 7i. Namely such operators are all candidates of annihilation- creation 
operators for 7i: 

k 

e iHt Ts{L ei L 62 ■ ■ ■ L ek )e~ iHt = Ts{L ei L 62 ■ ■ ■ L €k ) e iujmt , = m. (2.29) 

3=1 

Note that the order of L + and L~ is immaterial. 

The results and statements in this subsection are universal, i.e. they are valid for any 
root system A and for any choice of the Lax pair matrices, i.e. the choice of JZ. For deriving 
the explicit forms of the Heisenberg operator solutions, however, it is convenient to choose 
as 1Z the set of the least dimensions for each A. For the classical root systems, A, B and D, 
they are the set of vector weights^: 

A r _x : K = {ei,e 2 ,...,e r }, d = #n = r, (2.30) 

B r ,D r : TZ = {±ei,±e 2 ,...,±e r }, d = #K = 2r, (2.31) 

in which {ej} are the orthonormal basis of W, ej ■ e& = 5jk- We list in the Appendix the 
preferred choice of TZ for the exceptional and the non-crystallographic root systems. 



To be more precise, for B r they should be called the set of short roots. 
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3 Exact Heisenberg Solutions 



Now let us proceed to derive the explicit Heisenberg operator solutions for various 'sinusoidal 
coordinates' [TJ E] as many as the degree of freedom. For simplicity of presentation, let us 
fix the angular frequency of the harmonic confining potential as unity, u — 1 hereafter. We 
have already known one exact Heisenberg operator solution. In p] it was shown that the 
quadratic invariant 

r/W = J>J (3.1) 

is the simplest sinusoidal coordinate: 

[H,[H,v {1) ]} =4(?7«-c 1 W-c 2 ), Cl ,c 2 : const, (3.2) 

where c\ is given by rf^ = c\q 2 and c 2 = c\Eq. As shown in [TJ the exact Heisenberg solution 

e mt v (i) e -mt (3 _ 3) 

is easily evaluated. It should be noted that any root system has a degree 2 Coxeter invariant 
2 G -Fa (Table I), which is proportional to r]^\ In fact, 77 W is a sinusoidal coordinate [TJ for 
a Hamiltonian 7i ge more general than the Calogero system; a general homogeneous potential 
of degree —2 within a confining harmonic potential |13j : 



n 



^5>? + ^ 2 ) + ^)' £fc£n?) = -2V(<?), ( 3 - 4 ) 



3 =1 J= l ^ 



which satisfies [W ge , [W ge , g 2 ]] = 4(g 2 — 7i gc ). 

For the Calogero system based on any root system A, we will derive the explicit forms 
of the Heisenberg operator solutions: 

j m T)^e- m \ j = l,...,r, (3.5) 

for a complete set of sinusoidal coordinates defined by 

{r^W 2 ),...,?^}, V ®£Ta{Qf')=Tr{Qf'), f,eF A . (3.6) 



In (13. 6p . the matrix Q fl 2 . 2 5 [) is diagonal therefore its total sum (Ts) is the same as the 
trace (Tr). Let us note that any Coxeter invariant polynomial in {qj} can be expressed as 
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a polynomial in {7/ 1 ), r]( 2 \ . . . ,r)( r >}. It is easy to verify oc g 2 for an root system. The 
higher sinusoidal coordinates for the classical root systems are: 

r r 

A r ^ : r^oc^gf 1 , fc = l,...,r-l; ^oc^g,, (3.7) 

r 

S n £> r : r/^cx^gf, fc = l,...,r, (3.8) 
3=1 



\\<li (3-9) 



except for rf r ' for Z) r which reads 

i=i 

See the remark at the end of the section. As shown in Table I, all the integers {fj} are even 
for B r and _D r , except for f r — r for odd r of D r . This is also related to the fact that 

Tr(g 2 ^) = ^gf+ 1 + ^(-g,f +1 =0, 
3=1 3=1 

for the Lax operators based on the vector weights (I2.3ip . 

The sinusoidal coordinates take various different forms for the exceptional and non- 
crystallographic root systems. Note that the overall normalisation of {77^} is immaterial. 

The derivation of exact Heisenberg operator solutions is quite elementary. Let us in- 
troduce a generating function of the total sum of the homogeneous polynomials in L + and 
IT: 

G^\s) = Ts((L+ + sL~) f i), s G C, fj G F A , (3.10) 

which is a polynomial in s of degree fj 

Si 

G W W = EW»«'. b M = h;-i- ( 3 - n ) 
1=0 

The coefficient &/ j; / J _2« is the total sum of a completely symmetric product consisting of fj — I 
times L + and I times L~ and it can be explicitly evaluated. As shown in (12.291) . we obtain 

< »< h,j, *<" L 2l)t - (3.12) 

Namely bf.-j.-21 is either an annihilation operator (I > fj/2) being the positive energy part 
or a creation operator (/ < fj/2) being the negative energy part or a conserved quantity 
(/ = fj/2), the constant part. The conserved quantities are in general not in involution 

l>f^.bj,:n. /()• j?k. (3.13) 
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On the other hand, for s = — 1 we obtain 



and 



L + - L~ = 2iQ 



h 

c?W(-i) = (h^tbCqA) = (2z)V } = £(-l)%;/ 3 -i>. 

1=0 



W -=(2i)"*X;(-l) , 6/i ! /i-a- ( 3 - 14 ) 



i=0 



Thus we arrive at the main result of the paper; the complete set of exact Heisenberg operator 
solutions for the 'sinusoidal coordinates' {77^'}: 

h 

n&(t) d ^e mt r]^e- mt = (2i)"^ £(-1)' e^'" 2 ^, j = l,2,...,r. (3.15) 

This clearly shows that r)^(t) is a superposition of various sinusoidal motions. The trivial 
fact that these 'sinusoidal coordinates' commute among themselves 

is translated into the commutation relations among the annihilation-creation operators 

£ H-Ji-^ h h;h-*m] = 0. (3.17) 



I+m:flxed 



Among the annihilation-creation operators belonging to r/^\ the two extreme ones corre- 
sponding to I = and / = fj have a special meaning. They consist of L + (L~) only 

b fi . fi = Ts((L + )^), 6,,. i: = Tb{{L-)*), b\ rJi = bj . j (3.18) 

and commute among themselves 

as is clear from (I3.17p . These special annihilation-creation operators have been known for 
some time EJ Q21 QI] . 
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Since the Z-th term in (13.15P is annihilated by d/dt — i(fj — 21), we obtain a linear 
differential equation with constant coefficients satisfied by rfj> (t): 

fj , 

n (^-«)^ w (*)=°- ( 3 - 2 °) 



step 2 



Equivalently this can be rewritten as 



ft 

\\ (ad(W) - l) ■ r] {j) = 0. (3.21) 



step 2 



Here ad(H) denotes a commutator ad(H)X = f [H, X] for any operator X. For even fj case, 
the factor d/dt — i(fj — 21) for I = fj/2 could be omitted as it annihilates the conserved 
quantity. Then we obtain a differential equation of order fj 

fj , 

II (j f - il )-(v {j) (t)-2-H ff>0 ) = 0, (3.22) 

step 2 

fj 

J] (ad(W) - /) • - 2-% ;0 ) = 0, (3.23) 

step 2 

instead of fj + 1 for odd fj case. In [TJ , for a wide class of solvable quantum systems with 
one degree of freedom, we discussed the 'closure relation' 

[h,[h, v }}= v r (h) + [h, v ]r 1 (h) + R- 1 (h) (3.24) 

with R t (ny = Ri(H). By introducing a±(H) as Rx(H) = a + (H) + a -(H) and R (H) = 
—a + (H)ct-(H), this closure relation is rewritten as 

(ad(H) + a+(H)) (ad(H) + a -(H)) (ij + ^(H)" 1 R-i(H)) = 0. (3.25) 

Eqs. (13.211) and (13.231) are multi-particle generalisation of this relation. 

When a sinusoidal coordinate rft' is multiplied to an eigenvector n of H, 

fi 

(2i) VV» = E ( U^/:./: 2/ 0n, (3.26) 

the Z-th term belongs to the eigenspace of H with the eigenvalue £ n + fj — 2Z. Thus these 
fj + 1 terms are all orthogonal to each other. This is a multi-variable generalisation of the 
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three-term recursion relation of the orthogonal polynomials of one variable [lj. Through 

a similarity transformation in terms of the ground state wavefunction (f>o(q) (12. lip , let us 
define 

L = O 1 o L o O , Z± = cf)~ l o L ± o 0o = L ± zQ, (3.27) 

bf r ,fj-2i = 0o 1 ° b ff,fj-2i ° ^o, j = l,...,r, Z = 0,l,...,/j. (3.28) 

Needless to say the identity 

0o 1 o Ts(A n ) o 0o = Ts((0 o 7 1 o A o o r) 

holds for any matrix A consisting of operators. Then one can present the corresponding 
results for the multi- variable orthogonal polynomials {P n (q)} (I2.18P constituting the eigen- 
vectors: 

h 

(2iy^P n (q) = YSrVth-J*-* Pn(q)- (3-29) 

1=0 

The operators L, L and bf .j -21 ar e closely related to the Dunkl operators [61 fT5]. 

The annihilation-creation operators provide an algebraic solution method of the Calogero 
systems. The entire Hilbert space is generated by the multiple application of creation oper- 
ators on the ground state wavefunction 0o((?) = e w ^: 

r 

II II & /-?^°, Ho^N, (3.30) 

j=l 0<!</ J 

l=fj (mod2) 

whereas 0o is destroyed by all the annihilation operators 

6 /j;Z 0o = O, -fj<l<0, Z = / J -(mod2). (3.31) 

Obviously the above states (13.301) are over-complete and the orthogonality of various eigen- 
vectors belonging to the same degenerate eigenspace is not guaranteed. For example, a 
complete basis of the Hilbert space is given by using bf j j j only [6] 

r 

Y[b n f '. f . -0o, Vti,eN. (3.32) 
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Remark on D-type Theory As noticed above, the D-type theory requires separate treat- 
ment due to the special f r = r. We need Lax pairs based on the spinor and anti-spinor 
weights: 

K^{±\e 1 ± 1 -e 2 ±.--±\e r }, d*#H = r- 1 . (3.33) 

Here the the number of — signs is even (odd) for the spinor (anti-spinor) weights. They 
both form 2 r_1 dimensional representations of the Lie-algebra D r , and these representations 
are called minimal. The minimal weights for the A, D and E root systems have played an 
important role in constructing simple Lax pair representations [I]. 

For -Dodd we simply use the Lax matrix L, ^ (L,-.) based on either the spinor (sp) or 
the anti-spinor (as) weights and proceed in the same way as above: 

S )= T ^P)+ S Sp)) r )> SGC ' fr = r > 
r 

= ^b( sp)r . r „2is l . (3.34) 

1=0 

Then we obtain the expansion of 77^ in terms of the annihilation-creation operators and the 
corresponding exact Heisenberg operator solution 



G U)(-V = (2*) r Ts(Q[ sp) ) = (2i)V r) = ^(-l) Z 6 (sp)n ^, (3.35) 

1=0 

r 

77 W = (2i) _r ^(-l)'6( sp ) r;r _2i oc qiq2'-q r , (3.36) 
1=0 

v (r) {t) ^ v ( r) e -mt = m -rj2(-l) l b (sp)r ,._ 2l e^-^. (3.37) 

1=0 

For Z^even, the situation is slightly more complicated, partly because of the existence of 
another sinusoidal coordinate r)( r / 2 \ which is a Coxeter invariant polynomial of {qj} of degree 
r, too. Here we prepare the Lax matrices based on both the spinor and the anti-spinor 
weights L,^, since 

r] (r) oc Ts(Q( sp) ) - Ts(Q( as) ) oc q x q 2 •■■q r . (3.38) 
It is quite elementary to verify ( 13.361) and ( 13.381) . Let us introduce two 'generating functions' 

GfiOO = Ts((Lj p) + sL^y) = X> (sp)r;r _ 2 ^, (3.39) 

1=0 

= Ts((L+ s) +sL^r) = $> (as) , ;r _ 2 ^. (3.40) 



1=0 
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Next we introduce the sinusoidal coordinate r/^ as their difference at s 



-1, 



vi r) «« (2i)-( G w (_i) - = Ts(g[ sp) ) - Ts(g[ as) ) 



r 



(2i)-^(-l)' (6 (sp)r; 



r;r— : 




(3.41) 



We obtain the corresponding exact Heisenberg operator solution 



r 



r,W(t) ^ (^-^(-lj'hp)^ - & Mr;r - 2! )e^ ( . (3.42) 



1=0 

This completes the derivation of the exact Heisenberg operator solutions for the D-type 
Calogero systems. 

4 Summary and Comments 

A complete set of exact Heisenberg operator solutions, as many as the degree of freedom, is 
constructed for the Calogero systems based on any root system, including the exceptional and 
non-crystallographic ones. Based on the complete set, one can write down the Heisenberg 
operator solution e lHt A e~ lHt for any operator A expressible as a polynomial in the sinusoidal 
coordinates {r)^}. This is the first demonstration of the exact solvability of multi-particle 
quantum mechanics in the Heisenberg picture. At the same time, these Heisenberg operator 
solutions provide the explicit forms of various annihilation- creation operators, as the positive 
and negative energy parts. Their commutation relations are, in general, quite involved. As in 
the simplest case of degree one quantum mechanics [TJ, these sinusoidal coordinates and their 
expansion into the annihilation-creation operators provide the explicit forms of the multi- 
variable generalisation of the three term recursion relations for the orthogonal polynomials 
constituting the mult i- variable eigenfunctions. The derivation of the Heisenberg operator 
solutions is a simple consequence of the universal Lax pair, which manifests the quantum 
integrability of Calogero systems based on any root system. 

Let us conclude this paper with a few comments on possible future directions of the 
present research. For better understanding of multi-particle quantum mechanics in general, it 
is desirable to enlarge the list of exact Heisenberg operator solutions. The obvious candidates 
are: the Sutherland systems [7] having trigonometric potentials, various super-symmetric 
generalisations of the Calogero-Sutherland systems [3 [11], the Ruijsenaars-Schneider-van- 
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Diejen systems [16J which are 'discrete' counterparts of the Calogero-Sutherland systems, 
and the (affine) Toda molecules. 

The newly found annihilation-creation operators suggest an interesting possibility of in- 
troducing multi-particle coherent states as common eigenstates of certain annihilation oper- 
ators. It is a good challenge to construct explicit examples of such multi-particle coherent 
states having mathematically elegant structure and/or practical use. 

A completely integrable system, including the Calogero-Sutherland-Moser systems, has 
the so-called hierarchy structure. It is characterised by the existence of mutually involutive 
conserved quantities 

n 1 ,n 2 ,...,n r , [n j ,H k ] = o, j,k = i,...,r, (4.i) 

which could be adopted as independent Hamiltonians generating different but compatible 
time-flows; t±, t 2r ■ ■ , t r , as many as the degree of freedom. It is a good challenge to construct 
common Heisenberg operator solutions to all the flows of the hierarchy 

e <£;=i«itj fj(k) e -iJ2U H ^, k = l,...,r. (4.2) 
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Appendix: The Preferred Choice of 7Z 

Here we list, for the exceptional and the non-crystallographic root systems, the set 1Z to be 
used for the explicit evaluation of the Heisenberg operator solutions. They are of the lowest 
dimensionality. 

1. E 6 : The weights of 27 (or 27) dimensional representation of the Lie algebra. They are 
minimal representations. 

2. E 7 : The weights of 56 dimensional representation of the Lie algebra. This is a minimal 
representation. 

3. E$: The set consisting of all 240 roots. 
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4. Either of the set consisting of all 24 long roots or 24 short roots. 

5. Either of the set consisting of all 6 long roots or 6 short roots. 

6. ^{m)' The set consisting of the vertices R m of the regular m-gon 

R m = [ Vj = (cos(2A;7r/m + to), sin(2fc7r/m + t )) G M 2 | A; = 1, . . . , m}, (A.l) 
in which t = [ti jlm) for m even (odd). 

7. H%. The set consisting of all 30 roots. 

8. if 4 : The set consisting of all 120 roots. 
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